Superfluid shells for trapped fermions with mass and population imbalance 
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We map out the phase diagram of strongly interacting fermions in a potential trap with mass 
and population imbalance between the two spin components. As a unique feature distinctively 
different from the equal-mass case, we show that the superfluid here forms a shell structure which is 
not simply connected in space. Different types of normal states occupy the trap regions inside and 
outside this superfluid shell. We calculate the atomic density profiles, which provide an experimental 
signature for the superfluid shell structure. 
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The recent experimental advance on superfluidity in 
polarized ultracold Fermi gas has raised strong inter- 
est in studying thephase configuration of this system 

0,llllllllliIllallESE3 The experiments have 
suggested a phase separation picture with a superfluid 
core surrounded by a shell of normal gas 0, 0> 0- 
This picture has been confirmed also by a number of 
theoretical calculations of the atomic density profiles in 
the trap 0, 0, E3, E3 ■ As Feshbach resonances between 
different atomic species with unequal mass have been re- 
ported [12j, the next step is to consider the properties 
of strongly interacting Fermi gas with both mass and 
population imbalance between the two spin components. 
There have been a few recent theoretical works in this 
direction 0, Q] , with focus on the properties of a ho- 
mogeneous gas. 

In this work, we consider the properties of a trapped 
strongly interacting Fermi gas with mass and population 
imbalance. We map out its zero temperature phase di- 
agram in a harmonic trap (generally anisotropic) as a 
function of few universal parameters. Compared with 
the equal-mass case, the two-specie Fermi mixture ( 6 Li- 
40 K mixture for instance) shows a very different picture 
of phase separation: it supports a superfluid shell state 
in the intermediate trap region, with normal gases of dif- 
ferent characters filling the center and the edge of the 
trap. This unusual phase separation picture with non- 
monotonic superfluid order parameter in space only oc- 
curs for trapped fermions with unequal mass. We provide 
an intuitive explanation for the phenomenon, and show 
how to detect it by measuring the atomic density pro- 
files. This superfluid shell state is not simply connected 
in space, so it may support interesting vortex structure 
under rotation of the trap. 

Strongly interacting Fermi gas near a wide Feshbach 
resonance can be well described by the following single- 
channel Hamiltonian: 
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where e^ a = k 2 /(2m CT ) with m a denoting the mass of 
species-cr (or called spin a, with a —1,1 and h = I), 
V is the quantization volume, and aL is the fermionic 
creation operator for the kcr mode. The bare atom- 
atom interaction rate U is connected with the physical 
scattering length a s through the renormalization relation 
1/17 = 1/U P - (1/V) E k V(2ek) with U p = 4™ s /(2m r ) 
( m r = m^mi/(mi + mi) is the reduced mass). We take 
the local density approximation so that /if = fi T +h, fi[ = 
Mr - h, fi T = fi - V (r), where V (r) = J2 l=x , y ,z Pi T i/ 2 is 
the harmonic trapping potential, which, without loss of 
generality, has been assumed to be the same for the two 
components. The chemical potential fi at the trap center 
and the chemical potential imbalance h are determined 
from the total atom number N = + N± and the popu- 
lation imbalance p = (N-\ — N±) /N through the number 
equations below. 

Our calculation method is based on a simple exten- 
sion of the formalism described in Ref. [S| to the un- 
equal mass case. Under the mean-field approximation, 
the thermodynamical potential = — Tln[tr(e~ H / T )] (T 
is temperature) can be written as (lij 



n = -V|A| 2 /C/ + 5>k4-W-£ki] 
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7?k T (h + aekr) with 
k 2 /(4m r ), and 



E k = 

(m T - 



V( e kr - Mr) 2 + A 2 , ekr = K"/(4m r j, ana a 
mi) /(mi + mi). If one of the energies E±^ has one or 
two zero(s) in k-space, it signals the presence of the type- 
I or type-II breached-pair (Sarma) states 0, 0, H, E| 
(called the BP1 or BP2 states, respectively). The BP 
states represent a spatially homogeneous superfluid, but 
they differ from the conventional BCS states by a phase 
separation in the momentum space and by a topological 
change of the Fermi surface for the excess fermions. To 
determine the phase at different trap regions, we search 
for the global minimum of the thcrmo-potcntial Q with 
respect to the variational parameter A instead of using 
the gap equation dCl/dA = 0, as the latter may give 
unstable or metastable phases as remarked in Ref. 0]. 
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The \i and h are determined from the two number equa- 
tions N a — f d 3 rn ra integrated over the trap. The local 
atomic density n rai derived from the thermodynamical 
potential O as dfl/dfi a — —n ra V, has the expression 

n - = p Et M k/(^k, CT ) + vlf(~E k ,- a )}, (3) 

k 

where the parameters w k — (Ek + (ek r — Mr))/ (2£k), 
v£ = (Ek — (ek r — Mr))/ (2^k), and the Fermi distribution 
/(-B) = 1/ (l + e £ / T ). For convenience, we take — t=l 
and vice versa. The above mean-field formalism is also 
identical to the GqG diagram scheme if we interpret A at 
finite temperature as the total gap which includes con- 
tributions from both the order parameter and the pseu- 
dogap associated with the pair fluctuations || 0] . 

We calculate the phase diagram of the trapped 
fermions in terms of several dimensionless universal pa- 
rameters. For that purpose, the unit of energy is cho- 
sen to be the Fermi energy Ep at the center of the 
trap for N non-interacting fermions with an effective 
mass of 2m r and with equal population for the two com- 
ponents. Under the local density approximation, one 
finds Ep = {3Ny//3 x /3 y p z ) 1/3 /V^ from this defini- 
tion. The trapping potential V(r) has the dimensionless 
form V(r)/Ep = X^^fj where the normalized coordi- 
nates fi = ri/Ri and Ri = y/2Ep/ ' fii is the Thomas- 
Fermi radius along the «th direction. The momentum 
k and the temperature T are measured in the units of 
kp = y / 2(2m r )Ep and Tp = Ep/ks, respectively. The 
system properties then only depend on four dimension- 
less parameters kpa s , T/Tp, the population imbalance p, 
and the mass mismatch a. In the following calculation, 
we take a = 0.74 corresponding to the 40 K- 6 Li mixture, 
as the latter is the most likely experimental system for 
the two-specie Fermi gas. 

In Fig. 1, we map out the zero-temperature phase 
diagrams for trapped fermions as a function of the pop- 
ulation imbalance p at several characteristic interaction 
strengths kpa 8 . The system shows a rich picture of phase 
separation in the trap. First, on the BCS side of res- 
onance with (kpa s ) 1 = — 1 (Fig. 1(a)), even in the 
equal population case (p = 0), we cross four different 
phases from the trap center to the edge. The trap cen- 
ter is occupied by a normal mixture state with the heavy 
fermions ( 40 K) in excess (denoted as NMk), which is sur- 
rounded by a shell of BCS superfluid phase (denoted as 
SF). Further out, there is a shell of a normal mixture but 
now with the light fermions ( 6 Li) in excess (denoted as 
NMu). The trap edge is occupied by the single compo- 
nent normal gas of light fermions (denote as NPlO. The 
system behavior becomes significantly more c omp licated 
compared with the equal mass case 111 llll LLjj , where 
instead of several phases, there is only one superfluid 
phase over the whole trap in the corresponding configu- 
ration. Note that under the local density approximation, 
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FIG. 1: (color online) The zero temperature phase diagrams 
for trapped 6 Li- 40 K mixture near a Feshbach resonance. The 
phases include the BCS superfluid state (SF), the breached 
pair phase of type 1 (BP1), the normal mixture (NMk or 
NMLi, with 40 K or 6 Li in excess, respectively), and the nor- 
mal polarized (single-component) states (NPk or NPlO- The 
phases are shown versus the population imbalance p and the 
normalized trap radius r. 



as one moves out from the trap center, the chemical po- 
tential monotonically decreases, while the superfluid or- 
der parameter is apparently not a monotonic function in 
this case. The superfluid only occurs in an intermediate 
shell. We will give some explanations for this unusual 
phenomenon later. Continuing with the phase diagram, 
if we increase the population of the heavy fermions, the 
central NMk region grows, while the SF, NM^, and the 
NPtj phase regions shrink and finally all disappear at a 
critical population imbalance. After that, the trap edge 
is occupied by a single component normal gas of heavy 
fermions (denoted as NPk). If we increase the population 
of light fermions, the reverse happens. The central NMk 
region shrinks and finally disappear at a critical popu- 
lation imbalance, and the superfluid shell evolves into a 
superfluid core. 

On resonance (Fig. 1(b) with (kpa s ) 1 = 0), the 
superfluid phase region gets significantly larger. The 
normal mixture region NMLi at the intermediate shell 
completely disappears, which is quite different from the 
equal mass case where there is always such a mixed 
shell . Moving on to the BEC side of resonance with 
(kpa s ) = 0.2 (Fig. 1(c)), the type-I breached pair 
phase (BP I) appears at an intermediate shell between the 
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SF and NPu phases, but only when the light fermions 
are in excess. Compared with the equal mass case 0, 
the critical k F a s for the appearance of the BP1 state is 
significantly shifted towards the resonance point. An- 
other notable feature at this (fcj?a s ) _1 is that when the 
heavy fermions are in excess, the superfluid shell is not 
surrounded a normal gas anymore. All the normal com- 
ponents are pushed to the central core. Further on to 
the BEC side with (/Cira s ) _1 = 0.5 (Fig. 1(d)), the nor- 
mal mixture at the trap center finally disappears at all 
population imbalance, and we resume the picture of a su- 
perfluid core surrounded by a shell of a normal gas. The 
BP1 phase region at the intermediate shell grows as one 
expects, but again it only shows up when the majority is 
in the light fermions. 

A remarkable feature from the above phase diagrams 
is that the superfluid forms a shell structure in space, 
which separates different types of normal states at the 
trap center and the edge. This feature is qualitatively 
different from the equal mass case. Now we would like 
to understand in more detail how this feature shows up. 
We know that the phase is determined by the global min- 
imum of the thermo-potential O as a function of the gap 
A, under certain values of the chemical potentials \i v and 
h at the trap position r. As one moves out from the trap 
center, fi r monotonically decreases as fi r = fi — r 2 (in 
the unit of E F ) while h remains the same. In Fig. 2(a), 
we show ft as a function of A at several different val- 
ues of the normalized radius r (thus with different /i r ). 
The values of h and the central /j, are taken to be the 
typical ones for which there is a superfluid shell struc- 
ture in the phase diagram. The potential f2 typically has 
a double- well structure. At the trap center (the lowest 
curve with r = 0), the trivial well with A = is deeper 
which corresponds to a normal state. As one moves out, 
both wells are lifted, but with different speeds. At a lower 
critical value of r (which is 0.38 for the configuration in 
Fig. 2(a)), the two wells become equally deep. Above 
this value, the global minimum jumps to the nontrivial 
well with A ^ 0, which signals a first-order phase tran- 
sition to the superfluid state. As one moves further out, 
the nontrivial well approaches the trivial well; and at an 
upper critical value of r (which is 0.81 in Fig. 2(a)), 
the two wells merge, which signals a second-order phase 
transition from the superfluid to the normal state. Hence 
the potential f2 varies non- monotonically with /i r , which 
leads to the superfluid shell state only at the intermediate 
region. 

The above picture is established from the calculation 
of f2. We can also give an intuitive explanation for the 
superfluid shell state. Note that except for the deep BEC 
side with a very strong coupling, it is always more favor- 
able for the fermions to pair up when the mismatch of the 
Fermi surfaces of the two components becomes smaller. 
As one decreases the chemical potential /i r by moving 
out from the trap center, for non-interacting fermions 
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FIG. 2: (color online) (a) The thermodynamical potential 
Q shown as a function of the order parameter A (both in 
the unit of Ef) at different trap positions (different chemical 
potentials) characterized by the normalized radius r. The 
other parameters are p — 0.2, T = 0, and (fc^s) 1 = 0. 
The two solid curves bound a region corresponding to the 
SF phase, where the global minimum of Q is at a nonzero 
A. (b) and (c) Schematic illustration of the radius of the 
Fermi surfaces k F for the two components as a function of 
the normalized trap radius r. In (b), one can only have a 
superfluid core at the trap center, while in (c), one in general 
has a superfluid shell in the intermediate region. 



the radius of the Fermi surface in the momentum space 
k F (?) decreases as k F (r) = \/2m a (/v — r 2 ) (in the 
standard unit) for the component a. So, k F (?) for the 
heavy fermions decreases faster with increasing r. This 
qualitative statement should be true also for interacting 
fermions as interaction will not change the rough trend. 
Then, we can imagine two situations as depicted in Figs. 
2(b) and 2(c). If at the trap center, the Fermi surface of 
the heavy fermions ( 40 K) has a smaller radius k F (r = 0) 
(Fig. 2(b)), the mismatch of the two Fermi surfaces only 
grows as one increases r. Therefore the pairing super- 
fluid, if any, can only form a core at the center. On the 
other hand, if the Fermi surface of the heavy fermions 
has a larger radius k F (r = 0) (which is the case when 
40 K are in excess), the mismatch of the Fermi surfaces is 
minimized at an intermediate region (see Fig. 2(c)), so 
the superfluid only forms near that region and thus takes 
the shape of a spherical shell. This explains an impor- 
tant qualitative feature of the phase diagram in Fig. f . 
As one further moves to the BEC side, the Fermi surface 
mismatch (and thus the above mechanism) becomes less 
important, and finally becomes independent as to which 
component is in excess. The superfluid then always forms 
a core at the trap center where there is a larger atomic 
density (Fig. 1(d)). 

To detect the phase diagram of Fig. 1 in general and 
the superfluid shell state in particular, one can measure 
the atomic density profiles in the trap. The real-space 
density profiles for the polarized fermigas have been mea- 
sured in several experiments 0, 0, 1- Q: in particular, 
the most recent one has shown how to reconstruct the 
full density profile from the column integrated signal Q . 
We calculate the density profiles for several character- 
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FIG. 3: (color online) The atomic number densities n a (in 
the unit of np = fc|/(37r 2 )) shown versus the normalized trap 
radius r. The solid (dashed) curves are for the 40 K ( 6 Li) 
atoms, respectively. The inserts of (a) and (b) show the am- 
plified tails of the density profiles. The other parameters are: 
(a)T = 0, resonance, and p = 0.2, (b) T = O.ITf, resonance, 
and p = 0.2, (c) T = 0, (feo^ 1 = —I (BCS side), p = -0.4, 
(d) T = 0, (kpas)- 1 = 0.2 (BEC side), p = 0.3. 



In summary, we have mapped out the phase diagram 
of a strongly interacting fermion gas in a trap with both 
mass and population imbalance for the two spin com- 
ponents. As a remarkable feature, the superfluid forms 
a spherical shell in the case of mass mismatch. We at- 
tribute this phenomenon to a combined effect of the mass 
imbalance and the trapping potential. We also show that 
the superfluid shell should be clearly visible in experi- 
ments through explicit calculations of the atomic density 
profiles. The superfluid shell is not simply connected in 
space, so it should have nontrivial vortex properties un- 
der rotation of the trap, which is an interesting subject 
for future theoretical and experimental investigations. 

This work was supported by the NSF awards 
(0431476), the ARDA under ARO contracts, and the A. 
P. Sloan Fellowship. 



istic configurations of phase separation, and the results 
are shown in Fig. 3. Figure 3(a) is for the resonance 
case with a small population imbalance p — 0.2 and at 
zero temperature. The superfluid shell (where the densi- 
ties for the two components are equal) is clearly visible, 
which separates two normal regions. From the inside 
normal state to the superfluid, the heavy (light) fermion 
densities jump down (up), respectively. This jump is con- 
sistent with the first-order phase transition picture estab- 
lished from the thermodynamical potential f2 shown in 
Fig. 2(a). From the superfluid shell to the outside nor- 
mal regions, the atomic densities drop continuously (con- 
sistent with the second-order phase transition picture in 
Fig. 2(a)). There is a small region of the NM L i state, 
and outside is a tail for the NP L i region. In Fig. 3(b), we 
show the finite temperature density profiles (T = O.ITf) 
with otherwise the same parameters as in Fig. 3(a). The 
profiles get a bit more smooth (as one expects), but the 
jump from the inside normal to the superfluid shell is still 
clearly visible. Note that as pointed out in Ref. 0, the 
densities are not equal any more for the BCS superfluid 
state at finite T since quasiparticle excitations carry pop- 
ulation imbalance. Fig. 3(c) shows the density profiles 
on the BCS side. There are jumps in the density profiles 
from the superfluid shell to both the inside and the out- 
side regions of normal states. Fig. 3(d) shows the profiles 
on the BEC side, still with a superfluid shell, but there 
is no normal region outside the shell any more (no tails 
with unequal densities) . On the BEC side, there is also a 
region of the BP1 state at appropriate population imbal- 
ance. The BP1 state is hard to be directly seen from the 
real-space density profile, but can be unambiguously sig- 
naled with the detection of the momentum-space profile 
of the minority component fl5| . 
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